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We investigate mesoscopic fluctuations in the spin polarization generated by a static electric field 
and by Rashba spin-orbit interaction in a disordered 2D electron gas. In a diagrammatic approach 
we find that the out-of-plane polarization - while being zero for self-averaging systems - exhibits 
large sample-to-sample fluctuations which are shown to be well within experimental reach. We 
evaluate the disorder-averaged variance of the susceptibility and find its dependence on magnetic 
field, spin-orbit interaction, dephasing, and chemical potential difference. 
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A primary goal in semiconductor spin physics is the 
control of magnetic moments by electric fields [l], Q . One 
way to achieve this is to make use of the magnetoelec- 
tric effect (MEE) 0, ||, a spin polarized steady state 
which emerges from intrinsic 'magnetic' fields generated 
by spin-orbit interactions (SOI) and transport. While 
this MEE has been observed e.g. in n-InGaAs epilayers 
[1,0, 0] and hole gases H, the resultingnet polarization is 
below percent for electron systems [g, lZ| , and thus much 
smaller than what has been achieved by optical pump- 
ing [1, Moreover, in a standard two-dimensional 



electron gas (2DEG), with typical Rashba SOI [ll|, the 
MEE generates only in-plane spin polarization but no 
out-of-plane components Q. The latter would be desir- 
able, also since they can be detected more easily e.g. by 
optical means. 

However, these observations apply only to disordered 
phase-incoherent systems with self-averaging IT2 . On 
the other hand, it is well-known that phase-coherence in 
mesoscopic systems leads to new quantum effects such as 
conductance fluctuations or weak antilocalization, espe- 
cially due to intrinsic SOI 0, E3, El, H E3, 0] ■ 



Similarly, mesoscopic spin effects emerge for the MEE 
when the system becomes phase-coherent. Indeed, fo- 
cussing on 2DEGs with Rashba SOI, we will show here 
that the spin-electric susceptibility is subject to strong 
sample-to-sample fluctuations, and thus individual meso- 
scopic samples with phase coherence can exhibit a large 
net spin polarization. Quite remarkably, these strong 
fluctuations show up not only in the in-plane but also 
in the out-of-plane spin components. We find that 
these fluctuations considerably exceed the polarization 
obtained in self-averaging samples, and since the latter 
has been successfully measured by optical means [!, 0] , 
the spin fluctuations predicted here should be well within 
experimental reach. We will see that this strong enhance- 
ment is special for Rashba (or Dresselhaus[29j|) SOI, and, 
diagrammatically, it results from a spin vertex renormal- 
ization typical for such intrinsic SOI [3, O [2(| • 

Related effects studied before are local spin fluctua- 
tions in metallic conductors due to the extrinsic spin- 




FIG. 1: Dominant diagrams (for 1/pfI <C 1) leading to the 
variance b~^&X km given in Eq. see also Fig. [2] The upper 
diagram contains the Hikami boxes (HB) Vha, VR a (with solid 
arrows and two Diffusons (D)) and V^d,Vjid (with dashed 
arrows and two Cooperons (C)), resp. The lower diagram 
contains the HBs Vl!i,Vr& (solid arrows, two Diffusons) and 
Vlc,Vr c (dashed arrows, two Cooperons). 



orbit effect [2l| (as opposed to the intrinsic Rashba SOI 
22 1), density of states fluctuations of quantum corrals 
23j | , and fluctuations of spin currents in general nanos- 
tructures [24| and chaotic quantum dots 24 , [2B . While 
thermal spin fluctuations have been observed [261 ] , we are 
not aware of studies of mesoscopic spin fluctuations due 
to the MEE as described here. 

We consider a disordered mesoscopic square-shaped 
2DEG of size L 2 containing non-interacting electrons of 
mass m and charge e and described by the Hamiltonian 



H 



2m 



a{pia 2 - P2CT 1 ) + b • a + V. 
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Here, p = (pi,P2>0) is the in-plane momentum, a the 



Rashba SOI constant [ll|, 2b — g^s-B(cos ips, sin ipg, 0) 



an external in-plane magnetic field, and <r = (cr 1 , er 2 , er 3 ) 
the Pauli matrices (and cr° = 1). The disorder po- 
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tential V is due to static short-ranged impurities ran- 
domly distributed and characterized by the mean free 
path I = rpp/m, where r is the scattering time and pp 
the Fermi momentum. 

The spin polarization due the MEE is given in lin- 
ear response by (er 4 ) = X^Eji * — 1>2,3, where Ej is a 
static electric field applied along the j-direction and x y 
the (zero-frequency) spin-electric susceptibility (per unit 
area). Here, we focus on the mesoscopic fluctuations of 
X y due to disorder, described by the variance (Sx^) 2 , 
where 5x tJ = X y ~ X lJ an d where the overbar denotes 
disorder averaging. We start from the Kubo formula for 
X lJ expressed in terms of retarded/ advanced Green func- 
tions Gp at the Fermi energy Ep [lj| 



X ij (E f 



4tt 



Tra l (G| -G a Ef )v 3 (G 



E F ) i 



(2) 



where Tr — > J d 2 p/ (2ir) 2 trs denotes momentum integra- 
tion and spin trace, and Vj = i[H,Xj]/h is the SOF 
dependent velocity operator. In Eq. © we have used 
time reversal invariance [30| to make the symmetry x u = 



X Ji explicit [31 1 



The variance (Sx 1 -*) 2 is obtained as the impurity av- 
erage over the product of two susceptibilities given 
in Eq. ([2]). Extending the diagrammatic approach of 
0) EH 13 to include SOI and spin vertices, we obtain 
the diagrams shown in Fig. [1] which give the dominant 
contribution to the variance for 1/pFl <C 1. Explicitly, 
we find 



5 X »(E F + A)Sx km (E F ) = (^) 2 J 



d 2 q 



v l w:d^(A)d^(-a) 



(2tt) 2 

V^{D^(A)D^(A) 
+D^(-A)D«?(-A)} + V££Vfc { Cg M (A) C pK (q, A) 
+C7(-A)Cf (-A)} +V%V^C%"{-A)C«o{A) . 

(3) 

Here, A is the difference in chemical potentials or gate 
voltages of x lJ and x km \ an d D q and C q are the Diffuson 
and Cooperon matrices, resp., with 4x4 components 
/i, v = 0, 1, 2, 3 (see below). The Vlc's (Vr 's) are Hikami 



boxes (HBs) shown on the left (right) in Fig. [T] Since 
X H = x Ji ( se e Eq. ©), each product V£ p Vg" in Eq. © 
turns into a sum with 4 terms. These terms are obtained 
by exchanging spin and velocity vertices in the V's such 
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as e.g. Vj£ V^ a =V Lal u , u j v Ra yuj , 
(a 1 <-* «j) + (°" fc ^ w «j ' 0-1 ^ w j ' ) ■ Additionally, the 
vertices have to be dressed with non-crossing impurity 
line s ( 19l |20|). In contrast to conductance fluctuations 
T2I 14 1, such vertex corrections are crucial here as they 
give the dominant dependence on the SOI (see below). 

Let us now evaluate Eq. ([3]) by calculating first V L / R , 
given in Fig.[TJ and then C/D. From now on, we restrict 
ourselves to the diffusive regime 1/ L <C x — 2appT <C 1, 
which allows us to neglect the q-dependence in V L / R and 
to expand in x. Additionally, we may neglect b and A 
in Vt/fl. Indeed, we first note that Vi^o/H=o oc br and 
Va^q/Va=o oc At for small b and A, and second, that the 
suppression of C/D with increasing b and A sets in on a 



much smaller scale br 



and At « tf> (dephasing, 



see below). As a result, we can unify the calculation of 




FIG. 2: General Hikami box (HB): shaded area represents 
the correction to the 'empty' box V a by Vb,V c containing a 
single impurity line (dashed). The vertices are denoted by 
k 1 G {<r M , vi, V2}, and R/A stands for the averaged G 
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the 2x4x4 HBs in Fig. [T] They can be expressed by a 
linear relation 32j in terms of the general HB in Fig. [2 
defined by the 'empty' box Tr {G a k 1 G r k 2 G a k 3 G r k' 1 } 
and the associated HBs with a single impurity line. Here, 
k 1 G {<r^,vi,V2} denotes vertices in Figs. [IJ [2j and in 
Eq. ([3]), and G^^ b (p) is the impurity averaged Green 
function which depends on B-field and SOI [19]. 

Next, we evaluate the Diffuson D and Cooperon C in 
Eq. ©, given by D q = (l/2mr)(l - A+)" 1 and C q = 
(l/2mr)fc(l - X-^k, where k» l = tr{crW}/2. Ex- 
panding = Tr{^G R Fib (p)a-G A Fi±b (p - q)}/2mr 
in q = ql <C 1, /? = 2br <C 1, and x, we find 
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where </? is the polar angle of the momentum q and A = 1 — iAt. We will find below that the main contribu- 
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tion to Sx 2 comes from the 1/q 2 terms in D 00 and C 33 . 
To account for orbital dephasing we introduce a corre- 
sponding dephasing time in D and C by the standard 
replacement (see e.g. Eq.(3.15) in [13]) A — > A + i/r^. 

Although generalized to include spin vertices, the 
method presented so far involves the calculation of simi- 
lar diagrams as for the conductance fluctuations 12, 
An important difference, however, is the inclusion of the 
vertex corrections for spin and spin-dependent velocity 
which we discuss now. The vertex correction is an infi- 
nite sum of diagrams which consists of non-crossing im- 
purity lines connecting the retarded and advanced Green 
functions in either x % ^ or X km i resp. For the velocity ver- 
Pj/m, 



tex this leads to w,- 



i.e. the spin part of the 



velocity is cancelled in the dc limit [2C 
vertex this leads to the replacement a 1 



For the spin 
E'^cr^ where 

£ is diagonal given by £ = (1 — A + ) _1 at q — 0, with 
the relevant entries S 11 = £ 22 = 2/x 2 and £ 33 = l/x 2 . 
These expressions are valid in the regime l/L -C x and 
will be used here. For general x, L, the hnite size form of 
the vertex correction has to be taken into account, giving 
e.g. S 22 = (2/x 2 )(l-tanh(a;L/2/)/(a;L/2/)), which then 
renders Eq. ([7]) given below finite for x — > 0. 

In the regime max{(Z/L) 2 , t/t^} -C x 2 , we find for the 
in-plane (i = 1,2) and out-of-plane (i = 3) components 
of the variance 
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where vf = pp/m is the Fermi velocity, Sq = nl/L, and 
the sum over the q n = Sqn satisfies the mixed bound- 
ary conditions for the Diffuson [12j in a finite sample of 
square size L with two opposite sides attached to the 
leads. Here, s y in Eq. ([5]) are rational functions in 
Qn^yQn depending parametrically on x, b, A, and the 
orbital dephasing rate <f> — t/t^ = 2l 2 /L 2 ,. 

Evaluating s y for l/L,l/L$ <C x first for b = A = 0, 
and choosing the E-field along x (i.e. j = 1) we find 



ai (4xE") 2 
{q 2 + 2<t>) 2 ' 



(6) 



where a\ = a,2 = 
then obtain (Sx 



1\2 



1 and a 3 = 2. From Eqs. fSJ, ([S]) i 

(exT, m /2n 3 VF) 2 C2di, where a 



Stis n Sq i /[(nl+n 2 )Sq 2 + 2(j)}' 2 . To assess the magnitude 
of this result we compare it to the average of the in-plane 
susceptibility x 21 — ex/2-KVF For the out-of-plane 
component of the spin fluctuations this yields 



Zl\2 



2c 2 



(x 21 ) 2 



C2 T DP 

2vr 4 r 2 



(7) 



and similarly for the in-plane components. Thus, we see 
that the relative fluctuations grow with increasing top, 
where top = 2t/x 2 is the D'yakonov-Perel spin relax- 
ation time |28| . 



Remarkably, for negligible dephasing, i.e. <f> <C 1, one 
obtains C2 ~ 1.51, and the spin fluctuations (<5x ij ) 2 be- 
come independent of the sample size. Typical numbers 
for a GaAs 2DEG (that are consistent with the regime of 
validity for — 101) yield x = 0.1 ... 1 and thus result 

in large relative fluctuations, J (<5x u ) 2 /x 21 ~ 20 . . . 0.2. 
In other words, there exist specific impurity configura- 
tions and realistic system parameters that give rise to a 
large out-of-plane spin polarization in response to a static 
in-plane electric field. 
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FIG. 3: Example of a propagation path along which an ini- 
tially y-polarized spin (small arrow) rotates exclusively into 
the positive out-of-plane z-direction (denoted by circled dots). 
The spin directions are denoted by small arrows at the left of 
each segment and result from precession about the instanta- 
neous Rashba spin-orbit field. 

To gain physical insight into this result, we consider an 
electron with spin initially pointing along the y-axis, see 
Fig. [3] While the electron propagates coherently through 
the sample, the spin precesses about the intrinsic Rashba 
SOI field which is in-plane and perpendicular to the prop- 
agation direction. As a result of orbital phase coherence 
the electron propagates along a path that is preferred by 
constructive interference in the given disorder configura- 
tion. Fig.[3]shows an example of such a path and the spin 
directions associated with the propagation through each 
segment [33j|. Along this entire path the spin can only 
point up (+z-direction), but never down. Now, if initially 
the electrons were unpolarized, the net out-of-plane po- 
larization in this case would be cancelled by spins that are 
initially pointing along the negative y-direction. How- 
ever, due to the (in-plane) MEE, which itself is subject 
to strong fluctuations, e.g. due to conductance fluctua- 
tions, there is a finite in-plane polarization to begin with. 
The cancellation is therefore incomplete. These consider- 
ations make plausible that disorder configurations exist 
that give rise to strong out-of-plane spin polarizations. 

We next consider the effect of an in-plane magnetic 
field. For b = we see that the terms a) and d) in 
Eq. (|3|) contribute equally to the variance. However, for 
b ^ the 1 / g 2 -divergence is cut off in the Cooperon and 
we can approximate Eq. |5]) by making use of 
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1 + 1, 



32(6t) 2 
(g 2 + 2<f)x 12 



(8) 



where the first term in Eq. (J8]) results from the Diffuson 
contribution (the term containing Vhai^Ra m Eq. ([5])) 
which is not affected by (moderate) magnetic fields. 
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FIG. 4: The relative variance of the out-of-plane spin suscep- 

; — 2 

tibility (S\ 31 ) 2 / S\ 21 is shown as a function of the in-plane 
magnetic field 26r for x — 0.1, L<f, = 10CM, and sample sizes 
L = 1007rZ (solid line) and L = 2007r2 (dashed line), resp. 
The low field approximation, Eq. ([8j|, (curved dotted line) 
and half of the b = 0-value (straight dotted line) are shown 
for L — 7rZ100. Inset: Variance as a function of A for x — 0.1, 
L<f, = 100Z, and L = IOOttZ (solid line) and L = 200ttZ (dashed 
line). The suppression of <5x 31 with increasing 6 and A is 
described by Eqs. (0 and ®. 



Unlike the magnetic field, a difference in energies A 
(e.g. induced by gate voltages) leads to a suppression of 
all terms contributing to 5\ 2 - This is described by 



'A=0/ 



1 



4(Ar) 2 
(q 2 + 2<j>)< 



(9) 



which gives rise to a correlation scale for susceptibilities 
at different gate voltages. Indeed, according to Eq. ([9]), 
we can regard x(Ef) an d x{Ef + A) as uncorrected for 
A > max{0/r, (ttI/L) 2 /t}. 

In conclusion, we find strong mesoscopic fluctuations 
of the spin-electric susceptibility in a disordered 2DEG 
due to Rashba SOI, giving rise to a large out-of-plane po- 
larization. The predicted values and dependences on the 
SOI strength, B-field, dephasing rate, and Fermi energy 
are well within experimental reach. Such spin-dependent 
coherence effects, besides being of fundamental interest, 
might prove useful in spintronics applications aiming at 
the electrical control of spin polarization. 
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We note that Dresselhaus and Rashba SOI give the same 
results since the corresponding Hamiltonians are unitar- 
ily related 

The B-field breaks time-reversal invariance and leads to 
additional terms in Eq. {S). However, their contribution 
to (Sx'-') 2 is negligible for the small B-fields br ps ^Jx 2 (f> 
relevant for C, see Eq. ([8}. 

(G A G A ) drops out in Eq. J2J due to 

R/A _ ., f-fR/A-, 



the identity G^ A 



,G7 A }/h. 



E.g., Vg%[a%v m ] = (dk T ) Kl (dkk T yik» r V[a l ,ai,a r ,v m ], 
where k" 1 = tr{cr 2 cr M cr i }/2 and d{d) 
diag(l,-l,-l,-l,+l(-l),+l(-l)). 
Note that along each segment the electron undergoes 
many scatterings for x — l/lso *C 1- 



